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A SIMPLE MODEL FOR THE CLUSTERING OF SUBHALOS AS A FUNCTION OF MASS 



cn 

o 
o 

< 



> 

in 

m 

o 
m 
o 

Oh 
I 

o 



Michael R. Blanton^ 

Draft version February 2, 2008 

ABSTRACT 

Astrophysicists do not associate galaxies with the massive dark matter halos so easily detected in 
particle simulations and whose mass function can be predicted by the Press-Schechter formalism and its 
variants. Instead, they associate galaxies with smaller surviving subhalos whose contents do not fully 
mix into the larger dark matter halos they exist in. However, these subhalos are as of now difficult 
to find in simulations. I present a simple description of what the Press-Schechter formalism can say 
about the clustering of these subhalos, under the restrictive assumption that at fixed subhalo mass, the 
chance of survival of the subhalo is independent of its larger scale environment. This model predicts 
a weak dependence of clustering on mass below M « 10^^ A/0, in agreement with observations of a 
weak dependence of clustering on luminosity at luminosities less than L^. The model is in qualitative 
agreement with observations that at low luminosities galaxy clustering is a strong function of color while 
at high luminosities galaxy clustering is a strong function of luminosity. 

Subject headings: galaxies: clusters: general — galaxies: fundamental parameters — galaxies: statistics 
— large-scale structure of universe 



1. MOTIVATION 

The excursion set formalism of Press & Schechter (1974) 
(and the subsequent improvements of Bond et al. 1991, 
Lacey & Cole 1993, and Sheth et al. 2001) predicts the 
mass function of bound, virialized objects in the universe. 
The extensions of Mo & White (1996) and Sheth et al. 
(2001) predict the clustering relative to dark matter of 
these halos. However, the large-scale, sufficiently high res- 
olution iV-body calculations of Kravtsov & Klypin (1999) 
and other investigators show that each CDM dark matter 
halo can contain multiple subhalos, presumably ex-halos 
not fully mixed into the larger, surrounding halo. The 
halos themselves do not have properties similar to those 
of galaxies — most notably, they grow to be much more 
massive. However, it is plausible that galaxies reside in 
subhalos. 

Here I predict the clustering of subhalos as a function 
of redshift, mass, and scale, using a modification of the 
method of Mo & White (1996). I use the formation rate 
of halos to determine when to identify halos of a particular 
mass and then evaluate the clustering of those halos at the 
redshift of observation. Where cosmological parameters 
are necessary for these calculations I use those given by 
Bennett et al. (2003) in their Table 3. 

2. ESTIMATING THE CLUSTERING OF SUBHALOS OF A 
SPECIFIC MASS 

Consider a set of subhalos which when they reach mass 
M are identifiable galaxies that retain their unique identi- 
ties throughout history until z = 0, even after they merge 
into larger halos. Let us make one simplifying assumption: 
whether a subhalo of a given mass retains its identity dur- 
ing mergers into larger halos is independent of when it 
formed and its large scale environment. In other words, 
of all subhalos of mass M which ever existed, those which 
survive form a representative sample. This assumption is 
violated in the real universe and in simulations, but it pro- 
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vides an alternative to the standard treatment, in which 
each subhalo always loses it identity immediately after its 
formation. 

One can calculate the dependence on redshift of the for- 
mation rate of such subhalos from the extended Press- 
Schechter theory, the clustering at each redshift using the 
ansatz of Mo & White (1996), and integrate the results 
to the observed redshift Zobs using the continuity equa- 
tion (Nusser & Davis 1994; Fry 1996; Tegmark & Peebles 
1998). Notably, the same procedure cannot predict the 
abundance of subhalos as a function of mass, without fur- 
ther specification of the conditions of subhalo formation 
and survival (Sheth & Pitman 1997; Benson et al. 2001; 
Somerville et al. 2001). 

The Press-Schechter formulae are derived by assuming 
that the primordial overdensity field i5 is a Gaussian ran- 
dom field. The smoothed density field is distributed as a 
Gaussian with variance (T^(Af), where Al corresponds to 
the average mass contained within a smoothing length R. 
One associates each parcel of mass in the original distribu- 
tion to halos of mass Af at a particular redshift z if at that 
point in space, the smoothed density field on scale R lin- 
early extrapolated (using the growth function D{z)) to z 
is exactly some (5c, o- The conventional choice is (5c, o = l-^. 
We define 8c{z) — 5c.o/D{z). The fraction of mass in such 
a state is 



f{M,z)dM^ 
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X exp [-5l{z)/2<T^{M)] dM. (1) 

Lacey & Cole (1993) derive expressions for the proba- 
bility of merging between dark matter halos. The fraction 
of mass in halos of mass AI2 which just entered from halos 
of mass Ml is: 
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The formation rate of mass M halos is proportional to 
the fraction of mass in M halos which has entered recently 
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from all smaller halos. One can calculate this from the 
fraction entering M from each smaller mass halo, inte- 
grated over all possible smaller mass halos: 



To calculate the cross-correlation {5shS) I follow Mo & 
White (1996) and use the log-normal form for the distri- 
bution of the dark matter overdensity 5: 
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The divergence is irrelevant to the final result; below, I 
will discuss only the "shape" of the formation rate, which 
I express as 

S^(,)(-^)ey.p{^Sl/[2a'{M)]). (4) 



where x = ln(l + 6) + and ((5^) = exp(cr2) _ i. 

At linear scales the above equations can be reduced to 
a linear bias of the form 
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The Press-Schechter formalism has a simple prediction 
that the clustering of halos of a particular mass is inden- 
dent of their formation history. Thus the overdensity the 
newly born halos of mass M relative to their mean den- 
sity is the same as the overdensity of all halos of mass 
M relative to their mean, for which Mo & White (1996) 
obtain: 
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which has the notable property that it is always greater 
than unity. 

Note that I have not specified what overall number of 
subhalos end up at each mass! In fact, without further 
assumptions, one cannot do so. I have only assumed that 
the subhalos that survive at each mass have clustering 
properties representative of all the subhalos ever to have 
that mass. I ignore that this condition is almost certainly 
violated in the real universe. 

3. RESULTS 

First consider Figure 1 , showing the halo overdensity as 
a function of dark matter overdensity on scales of i? = 4 



_2cr2(Mo) 2a2(Mo' ' ' 

Assume a subhalo is in a mass overdensity d at the red- 
shift of observation Zobs) which determines the linearly ex- 
trapolated overdensity So- From this, one can find (5form(-M", 
{6h\S) at formation for subhalos of mass M formed at some 
redshift zj > Zobs which end up in mass overdensities of 5 
at redshift Zobs- Then: 
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To calculate quantities in the nonlinear regime, one must 
relate 5{) to (5, which I do using the formulae given by Mo & 
White (1996) for an Einstein-de Sitter universe (the solu- 
tion to the nonlinear spherical collapse problem is roughly 
independent of cosmology). 

The integrals above may be done explicitly to yield 

l + (4h(M,^obs)|5) = 

[l + ^(^obs)]exp(^g/[2a2(M)]) 

(72 (M) [l-a2(jy^^)/^2(^)]3/2 

x{A2exp(-(4-(5o)V[2A']) -S2exp(-; ' 
-F\/2^(SA2/t2),5o exp ($2) erfc(a;c)} 
where 

A2 = cr2(Af) -ct2(Mo) 
T2 = 4CT2(M)-3cr2(Mo) 
$2 = 252A2/[a2(Mo)T2] 
S2 = a2(Afo)[A2/T2] 

_ 4-(2a^(M)-a2(Mo))^o/T2 
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Fig. 1. — The relationship between overdensity of halos and over- 
density of dark matter from the treatment of Mo & White (1996) 

{on the left) and that from Equation 7 for subhalos (on the right). 
This figure has results for a smoothing scale R = 4 Mpc, for 
masses M = 10*, 10^°, and IO^^Mq, as labeled, and for redshifts 
2 = 0,1, 2, 3, 4, and 5, as labeled. 
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Fig. 2. — Top panels: Rate of formation dG/dz of subhalos of mass 
10^ Mq (left) and 10^^ Mq (^rij/iij in arbitrary units. Bottom panels: 
The solid line is the linear bias on large scales based on Equation 
10. The dotted line is the "bias at formation" from Equation 5. The 
dashed line is how you would extrapolate the results backwards from 
z = using only the continuity equation. The result of Equation 10 
is not far from assuming that all galaxies formed at near the median 
of the dG/dz distribution and were merely debiased from then on. 
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Fig. 3. — Bias at redshift z = as a function of mass for three 
different scales, as labeled. The solid lines are the predictions of 
Equation 7. The dashed liens are the predictions from Mo & White 
(1996) for halos identified at redshift z = 1.5 and observed at z = 0. 



/i^^ Mpc, for three different halo masses, at redshifts z = 
0. . . 5, as labeled. The left column is the results of Mo 
& White (1996) for halos identified at the same redshift 
as they are observed. The right column is the results of 
Equation 7. Particularly at low redshift, at fixed mass 
the relationship between (sub)halo overdensity and mass 
overdensity is much more linear and close to unbiased for 
the full set of subhalos than it is for the set of halos which 
just formed. This effect is commonly understood as re- 
sulting from the fact that the densest regions are usually 
fully contained in large mass halos by low redshift, and 
thus cannot form new subhalos. 

Consider now Figure 2. The top panels show the rate 
of formation dG/dz of subhalos of masses M — IQ^Mq 
on the left and lO^^M© on the right. The bottom panel 
shows as the solid line the results for large-scale, linear bias 
of Equation 10. The dotted line is the bias at formation 
of Equation 5. To demonstrate what is going on, I use 
the linear continuity equation to extrapolate backwards 
from the z — results of the solid line. This dashed line 
crosses the dotted line at around the median of the dG/dz 
distribution. The results for the bias at z = are similar 
what one would get by assuming that all the halos formed 
at the median dG/dz value. 

This description sounds similar to the usual method sug- 
gested by Mo & White (1996), to identify halos at one 
redshift and evaluate their clustering at another redshift. 
However, different mass halos form at different times. Con- 
sider the top panel of Figure 3, which shows as the solid 
lines the results of Equation 7 as a function of mass for 
smoothing scales of i? = 8 and 30 Mpc (as labeled), 
defining bias as 6 = (SshS) / {S5). For comparison, consider 
the dashed lines, which show the bias as a function of mass 
using the Mo & White (1996) results for halos identified 
at z — 1.5 and observed at z = 0. In general, for my re- 
sults the bias is less different than unity, because the bias 
near the median time of formation is typically closer to 
unity. Furthermore, the low mass halos which form very 
early naturally debias according to the continuity equa- 
tion. The bias is much less mass dependent under my 
assumptions than it is under the assumption that all halos 
stop merging at some fixed redshift. 

These results are extremely close to what one would get 
assuming all subhalos were formed at the median of the 
dG/dz distribution, determine by setting 



(52(z„ed) = 2a''{M)6l^ ln2 + (5,(zobs). 
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The dashed lines in the bottom panel of Figure 3 show the 
mass dependence using the standard Mo & White (1996) 
treatment but setting the redshift of identification at each 

mass to Zmed- 

Finally, although at low mass the bias is a very weak 
function of mass for the subhalos, recently assembled sub- 
halos are in lower density regions than older subhalos. Fig- 
ure 4 plots the bias of subhalos as a function of their mass 
and age, showing that at low masses age is important while 
at high masses mass is important. Accepting that galaxy 
luminosity is related to subhalo mass, and that the optical 
color of a galaxy relates to the time of subhalo assembly, 
this trend is qualitatively similar to the trend found in 
Hogg et al. (2003). Those authors show that for low lu- 
minosity galaxies, color is most closely related to galaxy 
overdensity, while for high luminosity galaxies, luminosity 
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is the most closely related to galaxy overdensity. 

However, the transformation from age of assembly and 
mass of a subhalo to color and luminosity of a galaxy is 
by no means simple or even necessarily one-to-one. In 
particular, the simplest interpretation of the hierarchical 
merging scenario of galaxy formation is that the most mas- 
sive galaxies have assembled the most recently; however, 
it is observationally the case that the most massive galax- 
ies have the oldest, reddest stellar populations. So if this 
scenario has any truth at all, it must be the case that age 
of assembly and color are not related at least for the most 
luminous galaxies, and so drawing firm conclusions for the 
low mass subhalos would be premature. 

Other results of Hogg et al. (2003) are more difficult to 
understand with this picture, such as the large overden- 
sities associated with red, low luminosity galaxies. This 
discrepancy may indicate that such galaxies preferentially 
survive in large mass clusters rather than in lower mass 
groups. 

4. DISCUSSION 

An observational consequence of the model presented 
here is that the clustering of galaxies should be a weak 
function of their mass, especially at low masses. Many 
papers have claimed that the average clustering is a weak 
function of luminosity at luminosities less than (Hamil- 
ton 1988; Giuricin et al. 2001; Norberg et al. 2002; Zehavi 
et al. 2002; Blanton et al. 2003; Hogg et al. 2003). Since dy- 
namical and lensing studies suggest that galaxies have 
total masses of order IO^^Mq (McKay et al. 2002), the ob- 
servational results are in agreement with Figure 3. In ad- 
dition, this model qualitatively explains the trend in Hogg 
et al. (2003) that at low luminosity, galaxy clustering is a 
strong function of color. 

These suggestive comparisons to observations inspire 
testing of the assumptions used here; with A^-body sim- 
ulations. Unfortunately, identifying subhalos in simula- 
tions is difficult. While many investigators are using high 
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resolution simulations and identifying subhalos, not much 
literature quotes directly subhalo bias with respect to mass 
overdensity. The results of Kravtsov & Klypin (1999) show 
that the relationship between their subhalo population and 
the mass overdensity is more linear and closer to 6 = 1 than 
the Mo & White (1996) prediction for halos, but that if 
they identify halos a,t z — Zobs + 1 and observe them at 
Zohs the formulae of Mo & White (1996) provide a rea- 
sonable fit. It would be interesting use similar simulations 
to predict clustering as a function of subhalo mass and to 
compare with the simple estimates I present here. 

Differences from this model in the observations or in the 
simulations may indicate a breakdown of the assumption 
that subhalo survival probability is independent of envi- 
ronment and may provide an insight into the nature of 
that process. 

A significant improvement of this model is to express 
the argument here in terms of the "halo model" of galaxy 
formation predicting the distribution of subhalos within 
each halo (see Berlind & Weinberg 2002 and references 
therein) . This requires an extension of the typical excur- 
sion set picture of Press-Schechter, and is implemented by 
Sheth (in preparation). 

I thank Ravi Sheth for helping me along with this project, 
David W. Hogg for useful discussions, and David H. Wein- 
berg for useful comments on the text. I am grateful for the 
hospitality of the Department of Physics and Astronomy 
at the State University of New York at Stony Brook, at 
which most of this work was done. 
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Fig. 4. Contour plot of the bias of subhalos b = {S^hS) / (SS) 
at redshift z = ou 8 h^^ Mpc smoothing scales as a function of 
mass and time of assembly. At low masses time of assembly is more 
strongly related to the large scale density field than it is at high 
masses. 



